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1. INTRODUCTION 

The Tropical Rain Measuring Mission (TRMM) 
conducted several intensive fieid validation 
campaigns for improved understanding of Tropical 
precipitation systems (Zipser et al., this 
proceeding). Two of the campaigns (TEFLUN in 
Florida and Texas, and LBA in Brazil) utilized: the 
NASA ER-2 high-altitude (20 km) remote sensing 
aircraft instrumented with the ER-2 Doppler Radar 
(EDOP), the University of North Dakota Citation 
microphysics aircraft, and the NCAR S-POL 
polarization radar. This paper focuses on EDOP- 
derived falispeeds and vertical velocities in the rain 
regions of two stratiform cases (5 September 1998 
along the east coast of Florida, and 17 February 
1999 in Amazonia in Brazil). These cases were 
sampled in situ microphysically by the Citation and 
reported elsewhere in this meeting; the main 
emphasis of this paper will be on the airborne radar, 
measurements and Inferences from them. 

2. BACKGROUND AND APPROACH 

EDOP is an X-band (9.6 GHz) nadir pointing 
Doppler radar. EDOP samples Doppler velocity, 
and reflectivity factor, Z, with a gate spacing of 37.5 
m (vertical resolution) and at 0.5 s intervals (100 m 
along-track sampling frequency). The antenna 
beamwidth is about 3°, a footprint of about 1.2 km 
at the surface. Calibration accuracy of EDOP is 
about 1 dB. EDOP provides similar measurements 
to ground-based vertically pointing radars except 
that extensive regions can be covered jn short 
times. 

Doppler velocities Vq from a nadir-pointing airborne 
radar represent the sum of the mean reflectivity- 
weighted hydrometeor motion, the vertical air 
motion, and the aircraft vertical motion. Atlas et al. 
(1972) derived relations between Doppler velocity, 
reflectivity, and rain rate assuming an exponential 
size distribution for rain. Uibrich (1994) expanded 
on this work assuming a Gamma size distribution. 
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Retrieval of information on raindrop size 
distributions with the above method requires that air 
and aircraft motions be removed from the observed 
Doppler velocities. 

The radar reflectivity factor Zin mm®m*^ is given by: 
Z=f^max/^(£>)£>6^£) 

where D (cm) is particle equivalent spherical 
diameter, N(D) in m'^cm'Ms number of raindrops 
per unit volume per unit size interval, D^is the 
median particle diameter. For a Gamma size 

distribution, N(D) = where p Is the 

shape parameter, and a=.A = 3.67 u * 10' 
o.3(9-}i)^ A/^( 6.4 X 10" exp(3.2p)) has been shown 

observationally and theoretically by Uibrich (1983). 
The integration limits and are assumed 
O-^oo ^ though several papers have dealt with 
truncated limits. The reflectivity-weighted Doppler 
velocity (ms’'') is given by: 

_ lgV(.D)N{D)D^dD 
\gN{D)D^dD 

where v(D) is the drop fallspeed. Furthermore, 
Uibrich (1994) has derived the reflectivity- weighted 
Doppler velocity v^as a function of p by assuming 
the Gunn-Kinzer fallspeed relation at ground level 
(1013 mb), i.e., v(D) = 9.65 - 10.3 e'^^. Assuming 
Wa\ vertical air motions are neglible, Eq. 2 can be 
given as a function of p and Z: 
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This relation applies to ground level (1013 mb); 
at other altitudes is obtained by multiplying it by 
[Po/pf where p and p^, are the air density at the 
surface and measurement height, respectively. Eq. 
3 has a major advantage over power law relations 
used extensively in the literature: is asymptotic 

with increasing reflectivity rather than increasing 
monotonically with Z. 


Figure 1 shows the dependence of v^on /musing 
Eq. 3; a commonly used power law relation is also 


plotted. At 40 dBZ, the variation of ^ from -2 to 6. 
results a fairly substantial 4 ms'^ variation in 
Figure 2 shows the Gamma size distribution 
corresponding to = 6, 8 ms*^ and =-2 to 4. The 
resulting range of generated curves are typical of 
reflectivities observed in EDOP data sets. Larger 
for a given ^ implies larger and a larger Z. 



0 20 40 60 

REFLECT;Vi’^ ^dBZ) 


Figure 1. v^- Z relation as a function of jj. 


3. EXAMPLES. 

If one has large sample of Vq and Z measurements 
from the same region, then g can be obtained by 
fitting Eq. 3. This requires the assumption that 
stratiform regions are fairly extensive. However, Eq. 
3 must be modified with an additive constant C 
since widespread vertical motions can be present in 
the stratiform region that will bias fits to Eq. 3. 
Smaller scale vertical motions can also be present 
in the stratiform but are Ignored in the following. 
Thus, the Vp- Z data from EDOP Is fitted with the 
following: 

xO.44 

— I v^ = C+v2(/i), (4) 

Po ) 

where is obtained from Eq. 3, C is a constant 
which consists of the mean mesoscale vertical 
motion, positive (negative) C indicates downdrafts 
(updrafts), and the altitude dependence of the 
measurements are removed by the density term on 
the left side. This fitting method provides a global 
estimate of p and C within some defined height 
interval, not a point-by-point estimate of the values. 
Besides the inhomogenlety 
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Figure 2. Gamma distribution for two values of 
assuming various values of the shape parameter p. 

of the size distribution over the analyzed region, a 
major source of error In this approach is the 
presence of small scale vertical air motions. In 
addition, small scale aircraft motions are not always 
fully removed from Vp resulting in an error in the 
calculation. 

Figures 3 and 4 show an examples of applying Eq. 
4 to stratiform regions in Florida during TEFLUN-B 
and Brazil during TRMM-LBA. The rain region 
between the surface and 400 m below the freezing 
level, is divided into approximately 500 m intervals 
and then fitted using a non-linear curve fitting 
routine. The vertical axis corresponds to observed 
Doppler velocities adjusted to the surface. The 
number of points used in Fig. 3 and Fig. 4 per layer 
are 6552 and 12000, respectively. It is noted that 
points above (below) the fitted curve are due to 
smaller scale downdrafts (updrafts) assuming 
aircraft motions have been removed. The scatter of 
the points around the fits is likely due to inadequate 
aircraft motion removal. Nevertheless, the fits are 
statistically significant, with an obvious increase in 
Vp with Z The fitted curves in Fig. 3 for stratiform in 
Florida indicate that granges from -1.4 at the top 
layer (2.74-3.19 km) to -0.4 at the bottom layer (0.3- 
0.75 km). Figure 4 indicates lower p values (-3) 
which suggest a narrower size distribution than for 
the Florida flight line in Fig. 3. The calculated C 
ranges from about 1. 5-2.0 ms ’ in the rain layers of 
Figs. 3 and 4 which suggests the presence of a 
mesoscale downdraft in both cases. Vertical 
reflectivity images for both cases suggest 
evaporation is possibly occurring for both cases 
which would account for the subsidence. 
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distribution of raindrops. This is in contrast to p>0 
reported in the literature. The radar measurernents 
and the approach presented will be explored further 
usinQ insitu microphysics data and pfound-based 
polarization radar data collected along with the 
EDOP measurements. 
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Figure 3. 

Observed measurements from stratiform 


on 5 September 1998 fitted with Eq. 4. 
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Figure 4. 

Same as Fig. 3 except for 21 February 

1999 fitted with Eq. 4. 





4. CONCLUSIONS. 

A "global" fitting procedure has been presented to 
obtain layer values of p in the Gamma distribution, 
and mean vertical air motions in stratiform rain 
regions. This approach facilitates a better 
understanding of the observational errors (aircraft 
motions) and the inhomogenieties In air motions 
and size distributions. The calculated pfor the two 
presented is less than 0, thus implying a narrow 




